function z = draw_frechet(E,epsilon,ndraws)

global ncity

rng(032878)
u = rand(ndraws,26);

for a=1:2
for d2=1:ncity+1
    z(:,d2,a) = (-(1/E(d2,a))*log(u(:,d2))).^(-1/epsilon);
end
end

